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The conduction of electrostatic fields from the outer magnetosphere to the ionosphere via the geo- 
magnetic field lines is studied in a quantitative fashion. It is assumed that any arbitrary pattern of 
turbulent fields in the magnetosphere can be analyzed into Fourier components, and the transfer of 
components whose wavelength lies between 0.1 km and 1000 km (as measured in the ionosphere) is 
examined. It is concluded that components of wavelength less than about 1 km will be severely atten- 
uated before reaching the F region, whereas components with wavelength greater than about 10 km 
will be only slightly attenuated. The application of these results to the theory of the origin of iono- 
spheric irregularities is discussed. 

The transfer of these electrostatic fields must necessarily be accompanied by the appearance of 
weak electrostatic fields along the direction of the magnetic field. The size of these fields is estimated, 
and their possible effect on the acceleration of electrons is discussed. 



1. Introduction 

The close relationship between the solar wind and the 
various types of geophysical activity peculiar to high 
latitudes, including the aurora, is now fairly well estab- 
lished. As a consequence, it seems reasonable to look 
for the basic cause of these effects in the dynamical 
interaction of the solar wind with the outer reaches of 
the earth's magnetosphere. Many of the consequences 
of this interaction have been outlined by Axford and 
Hines [1961] who suggested the existence of a large- 
scale convection of the material of the outer magneto- 
sphere driven by viscous forces at the interface between 
the magnetosphere and the domain of the solar wind. 
Because of the high electrical conductivity of the mag- 
netospheric plasma, the magnetic field is frozen to the 
material, and the convective flow pattern can be 
thought of in hydromagnetic terms as a steady motion 
of the lines of force throughout the magnetosphere 
and the major portion of the ionosphere. Alternatively, 
the motion can be described by saying that the motion 
of the material in the outer magnetosphere sets up an 
electrostatic field which exactly cancels the v X B 
dynamo field, thus allowing the material to flow freely. 
Because of the high conductivity along the direction 
of the magnetic field, the field lines are very nearly 
equipotentials of the electrostatic field, and the field 
is transferred down into the ionosphere, there driving 
the ionization in the direction of E X B. The exact 
equivalence of these two representations of the motions 
is a well-known result of the hydromagnetic theory 
of highly ionized plasmas, and is discussed in some 
detail by Axford and Hines [1961] and Hines [1963]. 

An important feature of the Axford-Hines theory is 
the suggestion that the outer magnetosphere is 



'turbulent', i.e., that perturbations of either v or B 
(or both) exist, and that much of the spatial fine- 
structure observed in high-latitude geophysical 
phenomena is due to the inward convection of this 
turbulence and its intensification by adiabatic com- 
pression of the tubes of force. There should be a 
corresponding turbulent component of the electro- 
static field appearing in the ionosphere, and this might 
lead to the formation of field-alined irregularities in 
ionization such as those giving rise to spread-/ 7 iono- 
sonde echoes, scintillation of radio-star signals, and 
auroral radar echoes. 

The mechanism is thus somewhat similar to that 
proposed by Dagg [1957] to account for the formation 
of F-region irregularities. He suggested that the 
upward transfer of turbulent electrostatic fields 
generated in the dynamo region of the ionosphere 
might lead to the formation of irregularities at much 
greater heights, and his proposal was examined in 
detail and criticized by Farley [1959, 1960, 1961], 
Dougherty [1959] and Spreiter and Briggs [1961]. 
The general conclusion was that turbulent fields of 
the appropriate scale-sizes would indeed be trans- 
ferred to the required altitudes, but that they would 
lead to the formation of irregularities only in the 
presence of an existing large-scale gradient of electron 
density, since the drift velocity in the F region must 
be very nearly divergence-free. Stated in another 
way, the random component of the electrostatic field 
will lead to a permutation of the tubes of force, but 
will lead to the appearance of irregularities only if the 
permuting tubes originally contained different amounts 
of ionization and are thus distinguishable from one 
another after permutation. The mechanism of Axford 
and Hines [1961] must operate under similar restric- 
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tions, and in discussing irregularities here we shall 
implicitly assume that the necessary large-scale 
gradients of ionization do in fact exist. The suggestion 
is then critically dependent on the efficiency with which 
small-scale electrostatic fields can be mapped down 
the field lines from the magnetosphere to the iono- 
sphere, and in this paper we attempt to examine this 
process quantitatively. One important consequence 
of the analysis to be described is the fact that the 
existence of turbulence in the magnetosphere must 
necessarily lead to the existence of weak electrostatic 
fields along the direction of the magnetic field in a 
steady state. Some of the consequences of the exist- 
ence of these fields will be discussed in a later section. 



2. Theory of Electrostatic-Field Mapping 

The geomagnetic field gives rise to a pronounced 
anisotropy in the electrical properties of the iono- 
spheric plasma, and the conductivity becomes a tensor 
quantity, i.e., an electric field vector E in a given 
direction gives rise to a current density vector j in a 
different direction. The anisotropy is most easily 
treated if we adopt a rectangular set of axes, one of 
which is along the local direction of the magnetic 
field vector, B. We shall assume here that B lies 
along the positive z-axis, and that the z-axis is directed 
vertically downward in the northern hemisphere. 
The current density is then given by (e.g., Baker and 
Martyn, 1953) 



j=^ 2 (E-B)B+|i{Bx(ExB)R^(BxE), (1) 

where o~o, en, 0*2 are respectively the longitudinal 
conductivity (along the direction of B), the Pedersen 
conductivity (in the plane of E and B and normal to 
B), and the Hall conductivity (normal to both E and B). 
Thus we have 



jx = (TiEx — <J<lEy 

jy = (TlEy + (TiEx 

j z =(ToE z . 



Since we are concerned with electrostatic 'fields, we 
can write 



E = — grad <p 



(3) 



where <p is a scalar potential, and in a steady state we 
assume 



divj = 0. 
Substituting (2) and (3) in (4), we find 



(4) 



dx 2 dy 2 



av, i d ( d<p\ n 



if we assume that the three conductivities vary only 
in the vertical direction. The Hall conductivity, cr 2 , 
does not appear in this equation, since our assumptions 
of a vertical magnetic field and a uniform ionosphere 
in the horizontal plane allow the Hall current to flow 
in uninterrupted divergence-free loops around the 
direction of the magnetic field, thus playing no part 
in determining the overall zero-divergence condition 
for the current flow. 

The problem is formally similar to that encountered 
in the calculation of electromagnetic fields in ground 
of inhomogeneous and anisotropic conductivity [e.g., 
Smythe, 1939; Wait and Conda, 1958] and was first 
treated for the case for the ionosphere by Farley [1959, 
I960]. We essentially follow the treatment of the 
latter author, who investigated the transfer of electro- 
static fields from the dynamo region of the ionosphere 
upward into the F region. In this approach the strong 
anisotropy in the conductivity is treated by a trans- 
formation of the coordinate system from (x, y, z) to a 
new system (£, 77, £), defined by 



£ = x 

rj = y 

d£ = (o-i/<ro) 1/2 rfz. 

Applying this transformation to (5), we find 



aV *V dV d£_a_ 

d? dy 2 •*? dt ^(ino-mj u, 

where 



(Tm = ((ToCT 



\l/2 



(6a) 
(6b) 
(6c) 



(7) 



(8) 



This equation is greatly simplified if we assume that 
(T m can be represented, at least over limited intervals 
of height, by the expression 



(2) (Tm = K exp (— pQ. 

Then (7) becomes 



ay ay ay d£ = 

de dr) 2 d£ 2 P d£ 



0. 



(9) 



(10) 



We now assume that the electrostatic potential is 
generated at a 'source level', £ , and has a sinu- 
soidal spatial variation with a wavelength A in both 
the £- and ^-directions. The sinusoidal variation can 
be regarded as one Fourier component of an arbi- 
trary distribution, and the assumption that the spatial 
wavelength is the same in both directions is introduced 
for convenience and does not affect the conclusions 
materially. Under these conditions, (10) is satisfied by 



<p = <p exp D4(£-£o)] exp 



2m 



(f + i?) 



]■ 



(id 
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where A is a root of the quadratic equation 

87T 2 

A*-pA-^ = 0, (12) 



i.e., 



2 



Pi 

4 



8tt 2 
" X 2 



(13) 



Let these two roots be denoted by 



.-§- 



p 2 877 2 ] 1 / 2 

.4 X 2 J 



i8 = 



[£+¥] 



1/2 



(14) 



The complete solution of (10) is then given by 
? = ¥>„[/» exp {«({-&)} + (> exp {/3(£ — Co)}] 



exp 



2m 



(£ + *?) 



(15) 



In order to make a close approximation to physical 
conditions, we divide the ionosphere below £ into 
n slabs, in each of which the conductivity a m is 
assumed to obey (9) and is continuous across the inter- 
faces between the slabs. Any actual profile of con- 
ductivity can then be used by adopting a sufficiently 
large value for n. The rath slab then lies between 
£ m _i and £ m , and the nth slab extends from £ n _i to 
infinity. Let the values of the coefficients in the mth 
slab be denoted by a„„ /3 m , P m , Qm- In the rath slab 
the potential then satisfies 

if = <p [P m exp {a m (£ - fm-i)} + Qm exp {frnit — im-l)}] 

\2<ni 



exp 



•(f + Tj) 



(16) 



Our boundary condition at the source level then 
leads to 



Pt + Qi = l. 



(17) 



At each of the (n— 1) interfaces between slabs the 
potential, <p, and its vertical derivative, d<p/d£, should 
be continuous. Thus at the interface between the 
rath and (ra-h 1) th slabs we have 

P m exp {a m (£ m — £m-\)} + Qm exp {j8 OT (£ m — £m-l)} 

= Pm + i + fe + i (18a) 

amPmexp {a^^- tm-\)} + PmQm e*P {j8m(£m — £m-l)} 
= a m+ l Pm+l +j3 m +i(?w» + i- (18b) 



The (n — 1) interfaces yield (2n — 2) equations; with the 
addition of (17) we have (2/i— 1) equations, and 2n 
unknowns (the values of P m and Q m in the /i slabs). 
The final equation is provided by adopting a boundary 
condition at infinity, the base of the nth slab. We 



assume that as £— > oo the potential vanishes. Farley 
[1959, 1960] also applied this boundary condition and 
was criticized by Spreiter and Briggs [1961] [see also 
Farley 1961] on the grounds that the conductivity of 
the ionosphere at great heights (£— > oo in Farley's 
model) is such that a vanishing potential leads to a 
finite current at infinity. The latter authors obtained 
solutions of (10) with the more realistic boundary 
condition of vanishing current at infinity. In the prob- 
lem considered by both Farley and Spreiter and Briggs, 
however, the direction of the positive £-axis was 
upward from the earth; in our case, the £-axis points 
downward, and the fact that the conductivity decreases 
rapidly to zero below the ionosphere leads to the dis- 
appearance of the current as well as the potential 
at infinity. Since a n must be negative and /3 n must be 
positive, this boundary condition leads to an equation 



Qn = 0. 



(19) 



The In equations given by (17), (18), and (19) now 
enable us to find the 2n unknown coefficients P m , 
Q m (m= 1, n). 

Having solved the equations, the 'transmission 
factor', Z, for the electrostatic potential can be found 
at any height from 

Z = <pl Vsource = Pm exp {a w (£ — £ m _i)} + Q m exp 

{/3m(£-£m-l)}. (20) 

This is equal to the transmission factor for the hori- 
zontal electrostatic field, since 

(!?)/(!f)' ooree = w/® 80 "" 6 " w * iouree - (21) 

The vertical component of the electrostatic field, 
along the direction of the magnetic field, is given by 



*■-(-)'"£■ 

\o- / d£ 



+ j8 w ()m exp {£„,(£-£„ 



/o-Al/2 

" ( — J <Po[<x m Pm exp {a,„(£ — £m-l)} 

m \ 27Ti 
i))]exp — 



(£ + !») 



We shall denote by L the magnitude of the ratio of the 
longitudinal component of the electrostatic field to 
the source potential at the top of the magnetic field 
line, i.e., 



L = l — I [a m P m exp {a m ((, — £m-\)} + fimQm 

• exp {)8«(£ -{«-i)}]. (23) 

It is worth noting that L is a maximum along the field 
line with the maximum source potential, whereas 
the horizontal field is a maximum where the gradient 
of the source potential is a maximum. Physically, 
this is due to the fact that the longitudinal field arises 
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from the necessity of supplying current along the mag- 
netic field to account for the horizontal variation of 
Pedersen current, and the latter is a maximum where 
the transverse field is varying most rapidly, not where 
it has its maximum value. 

3. Ionospheric Conductivities 

The conductivities cr and cri are functions of several 
quantities whose values can as yet only be estimated. 
The theory of the electrical conductivity of the ionos- 
phere below the peak of the F2 region has been 
reviewed by Chapman [1956], using as a basis the work 
of Chapman and Cowling [1939] and Cowling [1945]. 
For our purpose, we shall assume that the conductivi- 
ties are given to a sufficiently good approximation by 
simple mean-free-path theory. The longitudinal con- 
ductivity, o"o, is the same as the conductivity in the 
absence of a magnetic field, and is given by 

^ = ^4—+— 1, (24) 

in which the subscripts e and i refer to electrons and 
ions respectively, TV is the number density of ions or 
electrons, m and v are the respective masses and 
collision frequencies, and e is the electronic charge. 
We shall use MKSA units, in which e is measured in 
coulombs and the conductivities in mhos per meter. 
The electron collision frequency, v e , is assumed to 
contain contributions from both electron-neutral 
(v en ) and electron-ion (v e i) collisions, i.e., 



V e = V en + V e U 



(25) 



whereas v\ contains only a contribution from ion- 
neutral collisions, since the ions are almost unaffected 
by collisions with the much lighter electrons. 

Chapman [1956] derives the following effective 
values for the collision frequencies (converted to 
MKSA units) 

^ = 2.6 X lO-^/V+rc^r- 1 / 2 (26) 

v en = 5AXlO- 1 «n n T l l 2 (27) 

i/ ei = 4.2X 10- 6 {14.1 + log 10 (P/7V)}7Vr- 3 / 2 , (28) 

where W is the molecular weight of the ions, T is 
the temperature, and the subscript n refers to the 
neutral particles. 

The Pedersen conductivity, cri, is given by the mean- 
free-path theory as 



<n=Ne 2 



V e 



Vi 



m e (y\ + co 2 ) miv\ + cof)J ' 



(29) 



where (o e and m are respectively the electron- and 
ion-gyrofrequencies, i.e., 



0) e ,i z 



eB 

m,e,i 



(30) 




10* I0 3 10* I0 5 I0 6 

ELECTRON DENSITY (cm -3 ) 

FIGURE 1. Electron-density profiles adopted for daytime 
and nighttime conditions below 500 km altitude. 



Both conductivities depend on the electron density, 
both directly and through the electron-ion collision 
frequency as given by (28). For the purposes of this 
study we have adopted two profiles of electron density 
below 500 km, one appropriate to daytime and one to 
nighttime conditions. These are shown in figure 1. 

The conductivity relations given above are appli- 
cable only to regions in which there are many more 
neutral particles than ions and electrons (though 
v en is not necessarily large compared with v e i). At 
altitudes well above the peak of the F2 region this 
is no longer the case, and the conductivities approach 
those of a fully ionized plasma. Under these condi- 
tions, Chapman [1956] shows that the longitudinal 
conductivity is given by 



o-o = Ne 2 lm e Vei, 
and the Pedersen conductivity by 



0-1=0-0 



vL+<0* 



Substituting from (28) and (30), we find 
2.4X10VP/ 2 
/Me{14.1 + log 10 (P/AO} 



O"o = 



(31) 



(32) 



(33) 



ov 



4.2XlQ- 6 /7i e 

p/2 



{14.1+log 10 (F/A/)}, (34) 



since at great altitudes v\ { < a) 2 ,. 

In the present analysis we are not concerned with 
the values of the conductivities individually, but with 
their geometric mean, cr m , and with the depth-scaling 
factor (cri/cro) 1/2 . In the magnetosphere these are 



CFn 



N 



(35) 



and 
W^=-^^(f){14.1 + lo gl .(P/iV)} 



(=Veild)e). 



(36) 
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Figure 2. Height profiles of the geometric-mean conduc- 
tivity ((T m ) and the depth-scaling factor (cri/o" ) 1/2 for 
daytime and nighttime conditions. 

Both quantities are assumed to remain constant for all heights above 500 km. 



In both of these expressions, the height variation 
appears solely in the factor {N/B), if we ignore the slow 
variation caused by the logarithmic term in (36), 
and assume an isothermal magnetosphere. Data 
obtained from measurements of nose-whistlers, how- 
ever [Smith, 1961], have indicated that the ambient 
electron density in the magnetosphere is approximately 
proportional to the magnetic field strength, and the 
quantity (NIB) is thus approximately independent of 
height. In the calculations that follow, we have 
adopted a model for the magnetosphere in which (35) 
and (36) are applicable from the source height down to 
an altitude of 500 km; this is not strictly justifiable, 
but the likely departures from this height variation 
should not affect the conclusions materially. The 
value of p in this height range (see eq (9)) is then zero. 
If we assume that at an altitude where B is 10" 7 
Wb/m- 2 (100 gammas), T is 2000 °K, and TV is 10 9 
m -3 , we find that (o-i/o" ) ,/2 is about 4 X 10 -5 . In the 
analysis to be described, this value has been adopted 
for all altitudes in excess of 500 km. The resultant 
models for cr m and (o"i/or ) 1/2 are shown in figure 2. 

4. Conditions in the Magnetosphere 

We have shown above that conditions in the mag- 
netosphere (here defined rather arbitrarily as the 
region above 500 km) are relatively simple from the 
conductivity point of view, since both cr m and (o-i/o- ) 1/2 
can be regarded as constant with altitude to a first 
approximation. An additional complication is intro- 
duced in this region, however, by the rapid divergence 
of the field lines with increasing altitude. The model 
we have been discussing assumed that the magnetic 
field lines were parallel straight lines, a justifiable 
approximation within the ionosphere. The divergence 
of the field lines in the magnetosphere introduces two 
complications in the analysis. First, the spatial wave- 
length, X, increases continuously with increasing 
altitude, and second, the area through which a given 
current must flow along the magnetic field lines in- 
creases continuously with altitude, leading to a de- 
crease in the current density, j, with altitude for a 
constant total current. 



It can readily be shown that the separation between 
field lines in a dipole field is given by [e.g., Akasofu 
and Chapman, 1961]. 



s*« 



cos* 



^ 



(1 + 3 sin 2 V) 



1/2 



S^cos 3 ^, 



(37) 



(38) 



where ^ is the geomagnetic latitude of any point on a 
given field line. Here S x and S y are respectively the 
separations between points on a given pair of field 
lines in the ^-direction (normal to the field line and in 
the magnetic meridian plane) and the y-direction (nor- 
mal to both the field line and the meridian plane). For 
our purposes we require the separation as a function 
of distance along the field line, /, rather than of geo- 
magnetic latitude. The geometry of the dipole field 
leads to the following relation between the differentials 
of/ and ^: 



dl = -r e cos ^(1+3 sin 2 ^'W, 



(39) 



where r e is the radial distance of the field line in the 
equatorial plane of the dipole. Integration of this 
equation leads to 



/ = 



2V3 



[V3sin Wl + 3sin 2 ^o) 1/2 



V3 sin ^(1 + 3 sin 2 ^) 1 ' 2 



+ sinh- 1 ( V3 sin %) - shah" 1 ( Vs sin ^)] , (40) 

where / is the distance measured along the field line 
from a point with geomagnetic latitude ^ to a point 
with geomagnetic latitude ^. In what follows we 
shall measure / from the surface of the earth at a 
geomagnetic latitude of 65°, corresponding roughly 
to the location of the auroral zone. / is then given by 

/(km)- 10,270[4.152- {1.732 sin ^(1 + 3 sin 2 V) 1 ! 2 

+ sinh-Hl.732 sin ^)}] (41) 

and is illustrated in figure 3. The separation between 
neighboring field lines unit distance apart on the earth 
in the x- and y-directions is given by 






A sinusoidal electrostatic-field component of spa- 
tial wavelength X in the ionosphere now has spatial 
wavelengths Xj. and k y in the x- and y-directions in 
the magnetosphere, where k x , y = kS x ,y. The adopted 
invariance of <r m with altitude implies that p = in 
(13) and (14), and the term 87r 2 /X 2 in the same equations 
must be replaced by 477 2 (X~ 2 + X~ 2 ). 



24.58 cos 3 ¥ 
(1 + 3 sin 2 ^)" 2 


(42) 


13.25 cos 3 ^. 


(43) 
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FIGURE 3. Distance (1) measured from the surface of the 
earth along the dipole field line originating at geomag- 
netic latitude 65°. 



For the purposes of computation, the magneto- 
sphere is divided into a series of slabs, in each of which 
k x and k y are assumed to be constant, with discontinu- 
ous changes at the interfaces. Each slab now has 
unique values for a and (3 in (14), and Z and L can be 
found by matching boundary conditions across the 
interfaces. This model is open to question, since a 
correct matching across the slab interfaces can be 
achieved only at a limited number of points, and not 
over the entire x — y plane. Proper matching could 
be achieved only if the number of slabs used were 
infinite, in which case the variation of wavelength 
with height would be continuous. In practice a model 
consisting of 22 slabs was used, test computations 
having shown that the results differed negligibly from 
those of models using as many as 150 slabs, and thus 
by reasonable extrapolation from those of a quasi- 
continuous model. 

In the case of the ionosphere, the boundary condi- 
tions adopted were continuity of the potential, <p, and 
of the longitudinal electrostatic field, —(d<p/dQ. 

In the case of the magnetosphere, however, the 
decrease of current density with height introduced 
by the divergence of the field lines makes this latter 
boundary condition physically unsatisfactory since 
it would lead to a discontinuity in longitudinal cur- 
rent flow across the interfaces. A more desirable 
condition is the continuity of total current, i.e., of the 
quantity — SxSy(d(p/dQ, and the results to be de- 
scribed have been obtained using this condition to- 
gether with the continuity of the potential in the mag- 
netosphere. The condition expressed in (18b) is 
thus replaced by 



S m [a m P m exp {<*,„(£ — £ m -i)} 

+ p m Q m exp {j8 w (£-£ w -i)}] 
— S m +i[a: m +iP m +i + /3m+i(?m+iL (44) 

where the S's are the appropriate values of the prod- 
uct S X Sy. 



5. Method of Calculation 

The numerical calculations consist of the following 
steps: 

(1) The altitude range of interest is divided up into 
a large number of slabs. In those slabs lying above 
500 km altitude the conductivity cr m is assumed to be 
constant, but the spatial wavelengths of the electro- 
static-field components increase with altitude. In 
those slabs below 500 km, the spatial wavelength is 
constant with altitude, but the conductivity cr m 
obeys (9). 

(2) The conductivities are computed at the inter- 
faces between the slabs in the ionosphere. 

(3) The reduced depth, £ m , of each interface below 
the source level is found by stepwise numerical inte- 
gration of (6c). 

(4) The value if p is found for each slab below 500 
km from the computed conductivities at the top and 
bottom of the slab. 

(5) The appropriate values of a and /3 are found for 
each slab using (14) and assuming a value for the 
spatial wavelength, A, of the electrostatic-potential 
within the ionosphere. The wavelengths within the 
magnetosphere are found from the divergence of the 
dipole magnetic-field lines as described above. 

(6) The values of P and Q for each slab are found 
by solution of the simultaneous equations (17), (18), and 
(19) below 500 km, replacing (18b) by (44) above 
500 km. 

(7) The transmission factor, Z, and the longitudinal 
field-strength ratio, L, are found from (20) and (23). 

In practice, 67 slabs were used, 22 within the mag- 
netosphere (above 500 km), and 45 within the iono- 
sphere (below 500 km). The region from 500 km to 50 
km was divided into 44 slabs of 10 km thickness, and 
the final slab consisted of the region below 60 km. 
The COSPAR International Reference Atmosphere 
(CIRA, 1961) was adopted as a model for the neutral 
atmosphere below 500 km. The source level, at which 
the electrostatic potential originates, was taken as the 
equatorial plane, which lies 42,560 km above the sur- 
face of the earth, measured along the field line orig- 
inating at 65° geomagnetic latitude. The calculations 
were carried out on a digital computer. 



6. Results 

The results of the calculations are summarized in 
figures 4-11, and will be described briefly in this 
section. All values of the spatial wavelength, X, are 
to be understood as pertaining to the ionosphere. The 
actual wavelength at the source is considerably greater 
due to the divergence of the magnetic field fines. 

The potential transmission factor, Z, is shown in 
figures 4 and 5 for daytime and nighttime conditions 
respectively, and for a variety of values of the spatial 
wavelength, X. As would be expected intuitively, 
the downward transmission of electrostatic potentials 
is considerably less efficient by day than by night, due 
to the increased conductivity of the ionosphere, and 



832 




0.2 0.4 0.6 0.8 1.0 



Figure 4. The potential transmission factor (Z) as a 
function of altitude for daytime conditions. 

The labels on the curves refer to the spatial wavelength (A) in kilometers. 




Figure 5. The potential transmission factor (Z) as a 
function of altitude for nighttime conditions. 

The labels on the curves refer to the spatial wavelength (A) in kilometers. 



its consequently increased 'short-circuiting' action 
through the Pedersen currents. For both day and 
night conditions, however, it can be concluded that 
potential variations with wavelengths less than about 
1 km in the ionosphere are not likely to be transmitted 
from the outer magnetosphere, whereas those with 
wavelengths greater than about 10 km can be trans- 
mitted very efficiently. Thus potential variations 
over a global scale, such as those envisaged by Axford 
and Hines [1961] as arising from the convection of the 
magnetosphere, will appear in the ionosphere, and the 
motion of the ionospheric plasma will duplicate that of 
the magnetospheric plasma. The small-scale irregu- 
larities invoked by many authors to explain auroral 
radar reflections, however, are unlikely to arise in this 
way unless the corresponding potential variations in 
the magnetosphere are assumed to be of extremely 
large amplitude. Similarly, the irregularities thought 
to be responsible for radio-star 'fadeouts' [Little et al., 
1962; Moorcroft, 1963; Flood, 1963] are unlikely to 
arise in this way, since their spatial wavelength is 



thought to be considerably less than 1 km. The iono- 
spheric irregularities responsible for normal radio-star 
scintillation, however, have scale-sizes that are typi- 
cally of the order of a few kilometers, and may well 
arise from turbulent electrostatic fields transmitted 
downward from the outer magnetosphere, being attenu- 
ated in the process by a factor of 1 to 5. It is evident 
from figures 4 and 5 that the attenuation in the E region 
would be considerably greater than in the F region, 
so that this mechanism would not be expected to form 
such pronounced /^-region irregularities, especially 
during daytime. Observations (e.g., Yeh and Swenson 
[1964]) have shown that most, if not all, of the irregu- 
larities responsible for scintillation are indeed located 
in the F region. This point will be discussed further 
in the following section. 

The relative efficiency of transfer to the F region 
and the E region is illustrated in figures 6 and 7, 
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Figure 6. The potential transmission factor (Z) as 
function of A for daytime conditions. 
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The curves refer to altitudes of 500 km and 100 km. 
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Figure 7. The potential transmission factor (Z) as 
function of A. for nighttime conditions. 
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The curves refer to altitudes of 500 km and 100 km. 





which show Z as a function of X at 500 km and at 
100 km during daytime and nighttime conditions 
respectively. It is evident that electrostatic-field 
components with wavelengths lying between about 
1 km and 5 km can be transferred quite efficiently 
down to 500 km, but would not be expected to appear 
as low as 100 km. 

Figures 8 and 9 show the electrostatic-field compo- 
nent along the magnetic field which must necessarily 
be associated with the drop in potential illustrated 
in the previous diagrams. The parameter L is 
shown as a function of altitude for a variety of wave- 
lengths for daytime and nighttime conditions. It 
will be recalled that L is defined as the ratio of the 
longitudinal electrostatic field to the source-level 
potential at the top of the magnetic field line. In 
particular, if the amplitude of the sinusoidal potential 
variation at the source level is <po, then the maximum 
longitudinal electrostatic fields appear along the 
field line whose source potential is <po, and are given 
by —(d(p/dz) = L<po. 
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In the case of the smallest scale-sizes, these fields 
are confined to the magnetosphere and upper ionos- 
phere, since the electrostatic potential vanishes 
before reaching the F region. For the larger scale- 
sizes, L increases rather rapidly below about 200 km 
due to the increasing flow of Pedersen current and 
reaches a maximum well below the E region. The 
actual form of the height variation in this region will 
depend critically on the model adopted for the con- 
ductivity of the D region and the lower atmosphere; 
the model chosen here is not thought to be particu- 
larly accurate below 100 km, and detailed considera- 
tion of conditions at these altitudes will be deferred. 
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FIGURE 10. The longitudinal-field parameter (L) as a 
function of k for daytime conditions at altitudes of 
40,000 km, 500 km, and 100 km. 



Figure 8. The longitudinal-field parameter (L) as a 
function of altitude for daytime conditions. 

The labels on the curves refer to the spatial wavelength (A) in km. 
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FIGURE 9. The longitudinal-field parameter (L) as a 
function of altitude for nighttime conditions. 

The labels on the curves refer to the spatial wavelength (A.) in km. 



FIGURE 11. The longitudinal-field parameter (L) as a 
function of A. for nighttime conditions at altitudes of 
40,000 km, 500 km, and 100 km. 
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Since the motion of an electron along the direction 
of the magnetic field is unaffected by the presence 
of the magnetic field, these electrostatic fields may be 
expected to be important from the point of view of 
electron acceleration. In particular, those electrons 
in the high-energy tail of the thermal distribution func- 
tion may well take part in the well-known 'runaway' 
phenomenon of plasma physics [e.g., Rose and Clark, 
1961], gaining an energy comparable with the entire 
potential drop between the top of the magnetic 
field line and the altitude at which they make their 
first collision with a neutral particle. This suggestion 
will be discussed further, with some numerical esti- 
mates, in the following section. 

The variation of L with X is illustrated in figures 
10 and 11 for altitudes of 40,000 km (near the top of 
the magnetic field line), 500 km, and 100 km. At 
500 km Lisa maximum for a wavelength of 3 km during 
the day and 1.5 km during the night. At 100 km alti- 
tude the corresponding values of X are 5 km and 3 km. 
Thus the maximum ratios of longitudinal electrostatic 
field to source potential are attained in the ionos- 
phere for spatial wavelengths of a few kilometers, 
comparable with the scale-sizes of the radio-star- 
scintillation irregularities. For a given amplitude 
of source potential, figures 8 and 9 show that the maxi- 
mum total potential drop down to the F region is at- 
tained for values of X less than 1 km. However, as 
we shall see in the following section, the actual value 
of the source-potential amplitude is likely to be much 
smaller for these wavelengths than for the larger 
scale-sizes. 

7. Discussion 

It is of some interest to make a few numerical esti- 
mates of the sizes of the effects we have been discus- 
sing. As was mentioned in the initial discussion of 
the model, any arbitrary pattern of turbulent electro- 
static fields can be decomposed into Fourier compo- 
nents, and we have treated the components of different 
wavelength separately. The downward mapping 
process for an arbitrary pattern could be examined by 
recombining the Fourier components after multiplica- 
tion by the appropriate value of Z at any desired level 
in the ionosphere. In this discussion, however, we 
shall again confine ourselves to the consideration of 
individual sinusoidal components. 

A number of assumptions have been made in order 
to simplify the analysis; for instance, the magnetic 
field within the ionosphere has been assumed to be 
vertical, the conductivities have been assumed to be 
constant above 500 km, and horizontal variations in 
conductivity have not been permitted. Even with 
these assumptions, it is felt that the results have some 
quantitative significance. A more complete analysis 
would, of course, be highly desirable, but is outside our 
present scope. 

The large-scale magnetospheric convection system 
of Axford and Hines [1961] can be thought of as a 
Fourier component with wavelength of the order of 
1000 km in the ionosphere, and with a source-poten- 



tial amplitude, <£> , of about 10,000 V (derived from the 
observed ionospheric drift speeds of the order of 300 
m sec -1 and the known magnetic-field strength in the 
ionosphere). The results shown in figures 8 and 9 
then indicate that the maximum longitudinal electro- 
static field associated with this convection system 
would be of order 10~ 8 to 10~ 10 V/m -1 . Even if applied 
along the entire length of the magnetic field line, such 
potential gradients could only supply electrons with 
energies comparable with, or an order of magnitude 
greater than, their thermal energy. 

An estimate of the size of the source-potential ampli- 
tude associated with any component can be obtained 
from the magnitudes of ionospheric drift speeds. 
The electrostatic field, E, associated with a drift speed 
v is given approximately by vB. Also the maximum 
electrostatic field arising from a sinusoidal horizontal 
variation of potential of wavelength X is given by 
(27r/X)</?o, where <p is the source-potential amplitude. 
Equating these, we have 



</?o — 



kvB 

2tt* 



(45) 



The value of B in the ionosphere can be taken as 
3 X 10~ 5 Wb/m~ 2 ; appropriate values for v are largely 
a matter of guesswork, since v refers to the turbulent 
component of the drift speed, rather than its actual 
total magnitude. For our present purpose it might be 
reasonable to adopt values of about 100 m sec -1 and 
1 km sec -1 , the latter to be regarded as an upper limit, 
since it is comparable with the total drift speeds ob- 
served under fairly disturbed conditions. <p can 
then be taken as lying in the range 0.5 X to 5 X, where 
<£>o is measured in volts and X in kilometers. 

The longitudinal electrostatic field is given by the 
product Lipo, and is a maximum for wavelengths of 
the order of 10 km. Taking the upper limit for the 
source potential given above, we find (fo = 50 V at 
X = 10 km, and figures 8 and 9 then give a longitudi- 
nal field of about 5 X 10~ 7 V m _1 throughout the mag- 
netosphere. Thus the total potential drop along a 
magnetic field line from the equatorial plane to the 
top of the ionosphere is of the order of 20 V, a value 
in excess of the ionization potential of the atmospheric 
constituents. When we allow the possibility of an 
individual electron making several up-and-down trips 
in a region containing irregularities, it is apparent 
that these longitudinal electrostatic fields may well 
provide a source of energy for the ionization of the 
F region. 



8. Generation of the Fields 

We have considered the distribution of electro- 
static potentials necessary to provide a divergence- 
free flow of current transverse to the magnetic field 
in the ionosphere and along the magnetic field lines 
connecting the ionosphere to the magnetosphere, 
but we have not yet mentioned the remaining current 
component required to close the loop. Under certain 
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conditions (identical spatial variations of conductivity 
profile in both hemispheres), the completion of the 
current loop could presumably be provided by the 
ionosphere in the conjugate hemisphere without in- 
voking any contribution from transverse currents 
within the magnetosphere. In general, however, 
such a contribution must be supplied, and here we 
shall briefly discuss the nature of this current, which 
is intimately connected with the mechanism of genera- 
tion of the electrostatic fields. 

Any driving force acting on the magnetospheric 
plasma, such as the viscous force of the Axford- 
Hines theory, can be represented by a body force 
F acting on unit volume. Such a force acting on a 
highly ionized plasma causes the ions and electrons 
to drift in opposite directions, given respectively by 
the directions of the vectors FxB and — FxB. 
The resulting charge separation sets up an electro- 
static field E in the direction of — F X B, which in 
turn causes a drift in the direction of E X B, i.e., in 
the direction of the driving force F. If there were 
no ionospheric currents to consider, the field E 
would grow to a point at which the drift speed E/B 
equalled the solar wind speed. There would then be 
no net body force, and the 'polarizing' current would 
cease to flow, having achieved its purpose. In the 
presence of ionospheric currents, however, the polar- 
izing current must flow continuously to close the ionos- 
pheric-current loop, and enough resultant body-force 
must remain to generate this current. In other words, 
the magnetospheric plasma will accelerate only to 
a point at which the difference between its speed and 
that of the solar wind is sufficient that the resultant 
body force generates the required current. In this 
steady-state condition the current will exert a body 
force j X B on the medium, which will exactly balance 
the resultant driving force. The necessary current 
densities, j, in the magnetosphere are quite small. 
For example, an electrostatic field of 1.5 X 10 -2 V/m -1 
applied in the meridian plane throughout the daytime 
ionosphere will give an integrated Pedersen current 
of about 1.3 X 10" 2 A/m of longitude. If we assume 
that the driving force is applied over a length of some 
10,000 km at the outer edge of the field line, and allow 
a dilution factor of about 10 due to the divergence 
of the field lines, this current can be closed by a cur- 
rent density of about 1.3 X 10~ 10 A m -2 in the outer 
magnetosphere. If the electron density there is in 
the range 10 8 to 10 9 m~ 3 , this current density can be 
provided by a drift of ions and electrons in opposite 
directions with speeds of the order of 1 to 10 m sec -1 , 
which is extremely small in comparison with the drift 
speeds normal to E and B of about 10 km sec -1 re- 
quired by the observed ionospheric drift speeds of 
the order of 1 km sec -1 . Thus there appears to be 
no conceptual difficulty in visualizing the closing of 
the current loops by way of 'polarizing' currents trans- 
verse to the magnetic field in the outer magnetosphere. 
The full details of this process will not, however, be 
considered here. 



9. Conclusion 

We have attempted to examine quantitatively the 
ionospheric effects of electrostatic fields generated 
in the outer magnetosphere by motions of the plasma. 
It has been shown that the existence of small-scale 
'turbulence' in the magnetospheric motions may have 
far-reaching effects in the ionosphere, leading to the 
generation of irregularities in electron density, pri- 
marily in the F region, and possibly also to the accelera- 
tion of electrons to energies in excess of ionization 
potentials. 

As was mentioned in the introduction, however, the 
generation of F-region irregularities by this mech- 
anism requires the prior existence of gradients of 
electron density at the appropriate heights, and this 
point has not been discussed in the present paper. 
A further factor that influences the formation of ir- 
regularities is the lifetime of an individual electron; 
for example, an irregularity of one percent deviation 
from the background electron density can form in this 
way only if the average lifetime of an electron is suffi- 
ciently long that it can be transported by the random 
fields between points whose background electron 
density differs by one percent. The height of occur- 
rence of irregularities is thus likely to be determined 
by the most favorable combination of the various 
factors. 

The mechanisms involved appear to warrant fur- 
ther investigation, and in particular it would appear 
highly desirable to consider the effects of small-scale 
variations of conductivity in the horizontal direction, 
which have been ignored in the present discussion. 
During an auroral display, for example, the fine struc- 
ture observed in the visible light must be accompanied 
by corresponding horizontal fine structure in the 
ionospheric conductivity. Even if no turbulent struc- 
ture is present in the magnetospheric plasma motions, 
horizontal fine structure in the conductivity must give 
rise to horizontal variations in Pedersen-current flow, 
and corresponding longitudinal electrostatic fields 
associated with the large-scale convection of the mag- 
netosphere. A discussion of these effects will be 
deferred. It would also be desirable to investigate 
the conditions under which appreciable numbers of 
'runaway' electrons might be generated, and the possi- 
bility that these may contribute to the formation of 
the outer radiation belt. 

Yet another feature of the actual magnetosphere 
has been ignored, and merits further investigation. 
We have treated the geomagnetic field as though it 
terminated abruptly in the equatorial plane; in fact 
any given field line returns to the earth in the opposite 
hemisphere, and may there encounter an entirely 
different ionospheric conductivity profile. Thus any 
electrostatic field generated near the equatorial plane 
is faced with two 'circuits' in parallel, and the conse- 
quences appear to be well worth investigating. 
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